Two statistical-rigor gaps from review of scoring/stats.py. Both are analysis-layer and results-neutral for adherence — verified arms/cells/per_scenario/paired byte-identical before/after; only mcnemar gains additive keys and the Wald CI is clipped to its valid range. Full suite green (471 passed, 6 skipped).

1. Holm correction — promised but never implemented

spec/analysis-plan-amendment-1.md and paper/sections/method.tex both pre-register Holm correction for the H1 secondary cells, but no code computed an adjusted p-value — every McNemar p was rendered raw with no family-wise error control. This PR implements exactly what's pre-registered, no more:

• holm_adjust() — pure Holm-Bonferroni step-down (stdlib, deterministic).
• annotate_holm_family() — tags a crossover stats block per the frozen pre-registration:
  • rac-vs-naive_rag adherent at N∈{50,150} → p_value_holm + family="secondary" (Holm across whichever are present);
  • N=300 confirmatory → family="confirmatory", left uncorrected (the single pre-registered primary test);
  • everything else untouched (exploratory, reported raw and labelled).
• crossover.py routes its stats block through finalize_crossover_stats() at all three sites (single-seed, multi-seed, merged).
• Schema gains optional p_value_holm + family on the mcnemar object.
• Rendering: the report's paired table gains a Holm p / family column (secondary shows the Holm p, confirmatory reads "confirmatory (uncorr.)", exploratory shows "—"), a multiple-comparisons note, and the confirmatory sentence is now labelled "(uncorrected)" with a companion Holm summary for N∈{50,150}. The paper's stats_table.tex gets the matching column + footnote, so the table finally reflects method.tex's Holm claim.

Scope note: the family is exactly the pre-registered one — deviating to a broader post-hoc Holm would itself contradict the frozen plan. The H2 passed secondary (gitchameleon resolution path) is a separate family, an easy follow-on if wanted.

2. Wald risk-difference CI wasn't clamped

risk_difference(b, c, n) returned [diff ± 1.96·se] raw. A risk difference is bounded [−1, 1], but the Wald normal approximation can spill past it (risk_difference(9, 0, 10) → upper ≈ 1.09). Now clipped to the parameter space, so it never reports an impossible rate difference. wilson_ci was already clamped and the odds ratio is intentionally unbounded, so only this CI needed it. No schema change (the interval type is unbounded).

Testing

• holm_adjust on hand-computed vectors (m=1 identity, cap at 1, monotonicity, order preserved); annotate_holm_family correctly corrects the secondary family, leaves the confirmatory cell uncorrected, ignores non-confirmatory pairs, and handles partial N grids; a Holm-annotated crossover block validates against the schema.
• risk_difference clamp on out-of-range inputs (upper → 1.0, lower → −1.0); the existing in-range case still exact.
• Report + paper render tests for the Holm column, family labels, and the multiple-comparisons note.

Independent of PR #20 (roster Step 1); branched from main.